Superconductivity-Induced Distortions of Phonon Lineshapes in Niobium 
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Superconductivity-induced changes of phonon lineshapes in niobium have been re- investigated by 
high-resolution inelastic neutron scattering. We show that the changes go beyond a simple change 
in lifetime and frequency when the phonon frequency is close to the superconducting energy gap 2A. 
The observed lineshapes in elemental niobium are qualitatively similar to those found previously 
in borocarbide superconductors and agree very well with those predicted by the theory of Allen 
et al. [PRB 56, 5552 (1997)]. Our results indicate that the peculiar phonon lineshapes in the 
superconducting state predicted by the theory of Allen et al. are a general phenomenon and not 
restricted to a particular class of compounds. 

PACS numbers: 74.25.Kc,74.70.-b 



I. INTRODUCTION 



Superconductivity-induced changes of the frequency 
and the linewidth of phonons with strong coupling to 
the electronic system have been observed in conventional 
superconductors a long time ago^ 2 -. Similar effects were 
later found in high-T c cuprate superconductors as well^r— . 
When the phonon under investigation has a frequency 
smaller than the superconducting energy gap, there is 
a sudden drop in frequency as well as in linewidth on 
cooling through the superconducting transition temper- 
ature T c . This effect can be understood on the basis 
that the coupling of such a phonon to the electrons van- 
ishes in the superconducting state because phonons can 
no longer excite electron-hole quasi-particle pairs as soon 
as their energy is smaller than 2A(T). For phonons with 
frequencies above 2A an increase in linewidth was re- 
ported for T < T c (refs. [0|2|), which is due to the pile- up 
of electronic states just above 2A, i.e. the region in the 
electronic spectrum where the phonon couples to. 

The early publications on superconductivity-induced 
changes of phonon lineshapes concluded that these 
changes can be entirely described by changes of phonon 
frequencies and phonon lifetimes. In other words, it 
appeared that the intrinsic phonon lineshape remains 
Lorentzian below T c . However, marked deviations of 
the phonon lineshape from a Lorentzian were later ob- 
served in borocarbide superconductors, in particular in 
YNi 2 B 2 C (ref.i) and LuNi 2 B 2 C (ref.0). For a particu- 
lar acoustic phonon, the lineshape below T c did so little 
resemble a Lorentzian that the authors in Ref. [g thought 
to have observed a new type of excitations. These re- 
sults motivated theoretical work aiming at explaining the 
strange lineshapes observed in experiment. In one of the 
theoretical papers, i.e. that of Kee and Varma 8 , anoma- 
lous lineshapes were linked to a pole in the electronic 
polarizability in the superconducting state that appears 
near extremum vectors of the Fermi surface. This means 
that only phonons with wave vectors corresponding to a 
nesting vector of the Fermi surface are expected to show 
anomalous lineshapes in the superconducting state. The 



other theory put forward at that time, i.e. that by Allen 
and co-workers^, predicted lineshapes as observed in ex- 
periment quite generally, without any need for a strong 
nesting. 

The question which theory is more appropriate re- 
mained undecided for several years. Recently, we were 
able to show by studying phonons in YNi 2 B 2 C (ref. 
1 101 ) with various momentum transfers that the theory 
of Allen et ali^ gives a very good account of the experi- 
mental observations for all cases, which also means that 
nesting is not a prerequisite for the appearance of anoma- 
lous lineshapes in the superconducting state. Rather, the 
only relevant parameters are the ratios of the phonon fre- 
quency over the superconducting energy gap and of the 
linewidth in the normal state over the phonon frequency. 
Thus, similar effects as observed in borocarbide super- 
conductors should show up in other superconductors as 
well, and one might ask why they have not been seen so 
far. When comparing the phonon properties of various 
superconductors with those of YNi 2 B 2 C we found that 
the latter compound is particularly favorable for such a 
study. In most other compounds, the predicted effects 
are so small that they are easily overlooked. That is to 
say, the effects are nearly or completely washed out by 
the resolution of available neutron spectrometers. 

When searching for a superconductor other than the 
borocarbides to check the theory of Allen et al<2 we 
came upon elemental Niobium, which has a lower T c 
than YNi 2 B 2 C, i.e. 9.2 K instead of 15.2 K, but still a 
fairly strong electron-phonon coupling. The relatively 
strong electron-phonon coupling in Nb was the reason 
that it was chosen for an early inelastic neutron scatter- 
ing study of the effects of superconductivity on phonon 
frequencies and phonon lineshapes, i.e. that of Shapiro et 
ali£. The authors of this study did not mention to have 
observed any non-lorentzian lineshapes. It seems that 
the instrumental resolution in this study was insufficient 
to make any such deviation eye-catching. Very recently, 
superconductivity-induced changes of phonon linewidths 
have been investigated by a new technique, i.e resonant 
spin-echo neutron spectroscopy 11 . This technique is very 
powerful for detecting small intrinsic lineshapes, much 



smaller than can be detected by classical triple-axis neu- 
tron spectrometers. On the other hand, evaluation of the 
raw data relies heavily on the assumption that the intrin- 
sic lineshape is lorentzian in nature, and any deviation 
from a Lorentzian gives rise only to second order effects. 
Therefore, it is understandable that any deviations from 
a Lorentzian went unnoticed in this study. 

In this paper, we use again classical triple-axis neutron 
spectroscopy to investigate phonons in Nb above and be- 
low T c . We show that with improved resolution non- 
lorentzian lineshapes as predicted by theory are indeed 
observable. We further show that the theory of Allen 
et al. 9 gives a very good account of superconductivity- 
induced changes also for phonons with energies well be- 
low and above 2A, i.e. for phonons where the lineshape is 
not anomalous in the superconducting state, which fur- 
ther adds credence to this theory. 



II. EXPERIMENTAL 

The neutron scattering experiments were performed on 
the IT triple-axis spectrometer at the ORPHEE reactor 
at LLB, Saclay. Pyrolytic graphite was used both as 
monochromator and analyser. They were horizontally 
and vertically focusing, but tight collimation before and 
after the monochromator as well as before and after the 
analyser made the horizontal curvature ineffective. The 
tight collimation (25-20-20-30) was necessary to achieve 
a sufficiently good resolution in energy and momentum 
space. In principle, it would have been desirable not to 
make use of the vertical curvature as well in order to 
improve the momentum resolution in the direction verti- 
cal to the scattering plane. As will be discussed below, 
the relaxed vertical resolution led to asymmetric phonon 
lineshapes. However, tightening the vertical resolution 
would have led to exceedingly long counting times and 
therefore, we had to adapt our analysis to take the effects 
of the relaxed vertical resolution into account. Achiev- 
ing a good resolution required moreover to choose a low 
final energy, i.e. Ej = 5.0 meV. A graphite filter in the 
scattered beam was used to suppress higher orders. 

The sample was a single crystal of Niobium having the 
shape of a cylinder of 1.2 cm in diameter and 2 cm high 
with its long axis being parallel to the crystallographic 
(110) direction. Measurements were carried out in the 
100 — 011 scattering plane. The crystal was mounted in 
a standard orange cryostat allowing measurements down 
toT = 2K. 



III. RESULTS AND ANALYSIS 

Our measurements were restricted to transverse acous- 
tic phonons propagating along (100) for wave vectors be- 
tween q = 0.1 and 0.2 (in units 27r/a with a = 3.27 A). 
From the results published in Ref. [2j we knew that these 
phonons have a relatively large ratio of linewidth over fre- 
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FIG. 1: (a) Dispersion of the transverse acoustic phonon 
branch in the (100) direction. The line was calculated from 
the force constant model published in Ref. [l2|- Filled symbols 
depict our results for T = 10.5 K, open symbols are taken from 
fig. 9 in Ref. 2. (b) Phonon dispersion in the off-symmetry 
direction (0.16, %,x). The line was again calculated from the 
force constant model, whereas the open circles denote results 
of DFT calculations from Ref. [jj]. They were shifted up- 
wards by 0.9 meV in order to match the experimental value 
at (0.16,0,0). The arrow denotes the momentum resolution 
in the direction vertical to the scattering plane. 



quency and that the energies of these phonons would be 
comparable to the superconducting energy gap. Hence 
these phonons would be most favorable to search for 
superconductivity-induced changes of the phonon line- 
shape. The dispersion of these phonons is depicted in 
Fig. [T{a). These frequencies were extracted from energy 
scans like that depicted in Fig. [5J The asymmetry of 
the phonon line in Fig. [2ja) is due to the finite momen- 
tum resolution in conjunction with a steep increase of 
phonon frequencies when going away from the high sym- 
metry line (Fig. QIb)). The force constant model de- 
veloped by Reichardt and Pintschovius 12 by fitting the 
parameters to data published in the literature^ predicts 
a substantial frequency increase within the region in mo- 
mentum space sampled in our experiment (Fig. Hfb)). 
This prediction is supported by ab-initio calculations 
using density functional theory-is. We used the phe- 
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FIG. 2: (color online) Energy scans taken at Q = (0.14, 1, 1) 
at T = 2K (a) and 10.5 K (b). The thick (red) line in (a) 
represents a simulation based on the lattice dynamical model 
of Ref. [J2 and the resolution in Q and ui of the spectrome- 
ter. The resulting spectrum was fitted by 5 Gaussians being 
0.15 meV apart from each other (thin lines), (b) The thick 
(red) line was obtained by convoluting each Gaussian in (a) 
by the same Lorentzian and the same rigid energy shift to fit 
the normal state data. 



nomological model to simulate our energy scans based 
on the formalism of Cooper and Nathans^ for describ- 
ing the four-dimensional resolution function in q and en- 
ergy and found quantitative agreement with experiment 
(Fig. Ufa)). We note that the asymmetries related to 
resolution effects decrease with increasing q but are non- 
negligible up to the highest q values studied here, i.e. up 
to q = 0.2. 

In order to discriminate between asymmetries related 
to the finite momentum resolution and those related to 
superconductivity-induced distortions of the phonon line- 
shape, we adopted the following type of analysis: the 
phonon line predicted by the model was decomposed into 
five individual, equally spaced components having the 
linewidth expected for tight vertical resolution as illus- 
trated in Fig. [2] Subsequently, each component was con- 
voluted with the same Lorentzian to reproduce the line- 
shape in the normal state data. The Lorentzian linewidth 
obtained in this way (fig. [SJb)) was used as input for 



FIG. 3: Frequency (a) and linewidth (half-width at half- 
maximum) (b) of transverse acoustic phonons in the (100) 
direction observed above and far below T c . The solid line 
in (a) represents DFT results^ multiplied by a factor 1.5 to 
match the experimental data. The inset in panel (a) shows 
the observed energy shift AE on entering the superconduct- 
ing state, while the dashed, vertical line indicates the super- 
conducting energy gap 2A = 3.2 meV. The linewidth values 
were obtained by fitting the experimental spectra with a con- 
volution of a Lorentzian and the resolution function of the 
spectrometer as indicated in Fig. [2ja). Note that deviations 
from a Lorentzian in the superconducting state as discussed 
in the text were not taken into account in these fits. 



calculating the lineshape in the superconducting state 
within the framework of a theory published by Allen et 
alA The lineshape in the superconducting state pre- 
dicted by this theory was calculated for each component 
separately and the sum was used for comparison with 
experiment. For the phonon depicted in Fig. [5J theory 
predicts resolution-limited components somewhat down- 
shifted in energy. As will be discussed later, theory pre- 
dicts complex lineshapes for other phonons, i.e. those 
with energies very close the superconducting energy gap. 
For our analysis, we used the theory proposed by Allen 
ct al. 9 , because we had found that it works quite well 
for the effects seen in borocarbide superconductors^. 
It is based on the full quantum mechanical treatment 
of electron-phonon coupling where vibrational and elec- 
tronic excitations mix into hybrid modes. The theory 
contains three parameters, with two of them being en- 



200 




energy (meV) 



FIG. 4: (color online) Energy scans of a phonon with a fre- 
quency well above the superconducting energy gap 2A taken 
in the superconducting state (a) and in the normal state (b), 
respectively, at Q = (0.2, 1, 1). The red lines depict fits by 
convoluting the resolution function of the spectrometer with 
a Lorentzian. The blue line was calculated from the theory of 
Allen et al.-. 



tirely fixed by the phonon frequency and linewidth ob- 
served in the normal state. The third parameter is the 
superconducting energy gap 2A. The value of 2A was 
obtained by a comparison of the predicted to the ob- 
served phonon lineshape in the superconducting state. 
We note that the superconductivity-induced changes of 
the phonon spectra depend so sensitively on the gap value 
that the latter could be determined with a precision of 
about 1 %. 

In addition, we performed a simpler type of analysis to 
make contact with results of previous studies^ii: in this 
case, the phonon spectra were fit with a single Lorentzian 
convoluted with the resolution function. From these fits, 
we obtained the superconductivity-induced changes of 
phonon frequencies and phonon linewidths depicted in 

Fig. m 

The simple type of analysis just mentioned is fully 
adequate for phonons with energies well below or well 
above 2A: for phonons with energies well below 2A, the 
phonon lineshape becomes resolution limited below T c 
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FIG. 5: (color online) Energy scans taken at (a): Q = 
(0.16,1,1), (b): (0.17,1,1) and (c): (0.18,1,1) at tempera- 
tures above (red) and below T c (blue). The fit curves repre- 
sent a convolution of a Lorentzian with the resolution func- 
tion. The theoretical curves were obtained using the theory 
of Allen et al.— and a value of the superconducting gap as 
indicated by the dashed vertical line. The experimental back- 
ground is indicated by the thin horizontal lines. 



(Fig.EJa)), whereas for phonons with energies well above 
2A the intrinsic lineshape remains essentially Lorentzian 
(Fig. 2J. A close look at the shape of the phonon line 
depicted in Fig. HJa) reveals that there are indeed some 
discrepancies between the fit curve and the data, but they 
are small. Therefore, it would be hard to argue on the 
basis of such data that the lineshape in the superconduct- 
ing state deviates from a Lorentzian. Nevertheless, the 
curve calculated from the theory of Allen et al. 9 describes 
the data definitely better. 

We now turn to phonons with wavevectors q = 0.16 — 
0.18 where the lineshape in the superconducting state 



deviates very clearly from a Lorentzian. The energies of 
these phonons are quite close to the value of the super- 
conducting energy gap 2A = 3.2 meV. We note that we 
tried larger and smaller values of 2A as well but found 
the best agreement between calculation and experiment 
for 2A = (3.2 ± 0.02) meV, in agreement with the value 
reported in Ref. |2|. The lineshapes of these phonons are 
depicted in Fig. [SJ It is evident that the lineshape be- 
comes markedly more asymmetric in the superconduct- 
ing state, with a very steep intensity increase on the low 
energy side of the peak. Theory accounts both for the pe- 
culiar lineshape and the pronounced shift of the leading 
edge in a very satisfactory way. For a better understand- 
ing of the lineshapes predicted by theory, we show the 
theoretical curves before and after convolution with the 
experimental energy resolution in Fig. |S1 Although the 
fine structure of the theoretical curves is washed out by 
the finite resolution, our data lend strong support to the 
basic features predicted by the theory of Allen et al.— . 



IV. DISCUSSION 

It is obvious from Fig. Q] that the phonon frequen- 
cies obtained in our experiment match those reported in 
Ref. 13 extremely well. Further, we note that the phonon 
linewidths evaluated in a conventional way (Fig. |3|) agree 
with those reported in the two previous studies^ii very 
well. A close look reveals that there is some discrepany 
between our values and those reported by Aynajian et 
alf^i for the linewidths in the superconducting state 
around q = 0.16. Very probably, this discrepancy has a 
simple reason: that is to say, the measurements of Ayna- 
jian et al.— were done at a somewhat higher temperature 
than ours, i.e. 3.5 K instead of 1.8 K leading to somewhat 
larger linewidths. 

The aim of the present study was to show that the con- 
ventional analysis of the phonon lineshapes by fitting the 
data with a Lorentzian is not fully adequate if the phonon 
energy is very close to the superconducting energy gap 
2A. Our data obtained for wavevectors q = 0.16 — 0.18 
clearly show that this is the case, although the effects 
are less eye-catching than those previously observed in 
borocarbide superconductor s 6 ' 7 ' 10 . It is not that theory 
predicts much smaller effects for niobium than for the 
borocarbides, but they occur on a much-reduced energy 
scale and are therefore more washed out by the finite 
resolution of available neutron spectrometers. Neverthe- 
less, the very good agreement between our observations 
and the curves calculated from the theory of Allen et 
al. 9 confirm the predictions of this theory: if the phonon 
energy in the normal state is slightly above 2 A, but its 
low energy tail extends below 2 A, the low energy tail 
will condense into a very sharp peak right at 2 A on en- 
tering the superconducting state. According to Allen et 
al£, this sharp peak can be regarded as a mixed vibra- 
tional/superelectronic excitation. If the phonon energy 
in the normal state is right at 2A or even somewhat lower, 
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FIG. 6: (color online) The dashed (red) lines and the full 
(black) ones depict the theoretical curves before and after con- 
voluting them with the experimental energy resolution. The 
value of the superconducting energy gap used in the calcula- 
tions is 2 A = 3.2 meV. For the sake of comparison with fig. [S] 
the experimental background indicated by the horizontal dot- 
ted line has been added. The occurrence of two peaks in the 
top panel is an artifact related to the decomposition of the 
phonon line into a small number of individual components: a 
decomposition into a larger number of components would lead 
to densely spaced sharp peaks between 3 meV < E < 3.2 meV. 
However, the lineshape would barely change after convolution 
with the experimental energy resolution. 



the sharp peak is located somewhat below 2A and the 
high energy tail above 2A has less weight. If the phonon 
frequency in the normal state is considerably smaller than 
2 A, the high energy tail becomes vanishingly small, and 
the sharp peak corresponds to a normal phonon with in- 
finite lifetime. If the phonon energy in the normal state 
is considerably larger than 2A, there is still an intensity 
jump and a sharp peak right at 2A but the jump height 
is so small that it is very difficult if not impossible to see 
it in experiment. Our data for q = 0.20 shown in Fig. 0] 
do support the presence of such an intensity jump at 2A 
but only marginally so, the jump height being simply too 



small. 

We emphasize that the theory of Allen et al£ does not 
only explain the peculiar lineshapes observed for phonons 
with energies very close to 2A but also describes the 
superconductivity-induced frequency shifts and linewidth 
changes for phonons with smaller and larger frequencies 
in a quantitative way. We note that the agreement be- 
tween theory and experiment in Niobium is better than 
previously found in YM2B2C (Ref. 10). Possibly, the 
theory is better suited for Nb because it was based 
on a weak-coupling BCS approach and the electron- 
phonon coupling strength is indeed weaker in Nb than 
in YNi2B2C. Another reason might be the fact that the 
superconducting energy gap shows a marked anisotropy 
in YNi 2 B 2 C (Ref.Hoh but not in Nb, and that the theory 
of Allen et al£- was developed for the case of an isotropic 
gap. 



V. CONCLUSIONS 

We re-investigated the superconductivity-induced 
changes of the phonon lineshape in elemental Niobium 
using high-resolution inelastic neutron scattering. Our 



data confirm the results of previous studies but also show 
unambiguously that for phonons with energies very close 
to the superconducting energy gap, the effects of super- 
conductivity go beyond a simple change in lifetime. Our 
results strongly support the view of Allen et al£ that, 
in contrast to a competing theory-^, the peculiar phonon 
lineshapes observed in borocarbide superconductors are a 
general phenomenon and are not restricted to the case of 
a strong nesting effect. The effects found in Niobium are 
indeed quantitatively described by the theory published 
by P. B. Allen et al.— . The agreement between theory 
and experiment for Nb is better than for the borocar- 
bides. The theory does not only give a very good account 
of the peculiar lineshapes of phonons with energies close 
to the superconducting gap, but also of the changes in 
frequency and linewidth for phonons with significantly 
smaller or larger frequencies. The solid theoretical ba- 
sis of the superconductivity-induced changes of phonon 
lineshapes allows for a precise determination of the su- 
perconducting energy gap in a single phonon scan if the 
phonon energy is close to the gap energy. 

Work at Argonne National Laboratory was sup- 
ported by the U.S. Department of Energy, Basic Energy 
Sciences-Materials Sciences, under Contract No. DE- 
AC02-06CH11357. 
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